Quark-Gluon Plasma Fireball 



O 

o 
o 

(N 
^— > 
O 

O 

-xf : 

(N . 

: 

> : 

in . 
oo ■ 

o ■ 

: 

o . 
o ■ 
o ; 

^ : 

Oh! 

i 

a*: 



X 



Salah Hamieh, Jean Letessier and Johann Rafclski 
Department of Physics, University of Arizona, Tucson, AZ 85721 
and Laboratoire de Physique Theorique et Hautes Energies, Universite Paris 7, 2 place Jussieu, F-75251 Cedex 05. 

(Received: June 17, 2000) 

Lattice-QCD results provide an opportunity to model, and extrapolate to finite baryon density, the 
properties of the quark-gluon plasma (QGP). Upon fixing the scale of the thermal coupling constant 
and vacuum energy to the lattice data, the properties of resulting QGP equations of state (EoS) 
are developed. We show that the physical properties of the dense matter fireball formed in heavy 
ion collision experiments at CERN-SPS are well described by the QGP-EoS we presented. We also 
estimate the properties of the fireball formed in early stages of nuclear collision, and argue that 
QGP formation must be expected down to 40A GeV in central Pb-Pb interactions. 



PACS: 12.38.Mh, 12.40. Ee, 25.75.-q 



I. INTRODUCTION 

It is believed today that a new state of matter has 
been formed in relativistic nuclear collisions at CERN 
0. The existence of a novel non-nuclear high temper- 
ature phase of elementary hadron matter, consisting of 
deconfined quarks and gluons, arises from the current 
knowledge about quantum chromodynamics (QCD), the 
theory of strong interaction. At sufficiently high temper- 
ature, the strong interactions weaken (asymptotic free- 
dom), and thus we expect that hot and dense nuclear 
matter will behave akin to a free gas of quarks and glu- 
ons Q. At issue is today if the new phase of matter 
observed at CERN is indeed this so called quark-gluon 
plasma (QGP) phase. One of the ways to test, and pos- 
sibly falsify, the QGP hypothesis is to consider if the 
expected properties of the QGP indeed agree with ex- 
perimental data which have been at the center of the 
CERN announcement. 

The observational output of experiments we consider 
are particle abundances and particle spectra. We address 
here results on hadron and strange hadron production at 
the equivalent center of momentum (CM) reaction en- 
ergy Eqm = 8.6 GeV per baryon in Pb-Pb reactions. As 
seen in many experimental results [0 , which we will not 
restate here in further detail, in these high energy nu- 
clear collisions a localized dense and hot matter fireball 
is formed. In our earlier analysis of experimental results 
P,0 , we have obtained diverse physical properties of the 
source, such as energy per baryon content E/b, entropy 
per baryon content S/b, and last, not least, strangeness 
content per baryon s/b. These values are associated with 
temperature T/ and baryo-chemical potential /ib at which 
these hadron abundances are produced (chemical freeze 
out), along with other chemical properties and the col- 
lective velocity of the matter emitting these particles v c . 

Our main aim is thus to compare the properties of the 
QGP phase, modeled to agree with the lattice-QCD cal- 
culations to the properties of the fireball obtained 
in the study of hadronic particle abundances. We thus 
develop semi-phenomenological QGP equations of state 



(EoS) based on thermal and lattice-QCD results. These 
are then applied to explore, in a systematic fashion, prop- 
erties of the fireball formed in the Pb-Pb reactions. We 
show that the physical properties of the fireball extracted 
from particle production data are the same as obtained 
employing our QGP-EoS. 

Using our QGP-EoS, we also explore the initial ther- 
mal conditions reached in the collision both in, and out, 
of quark chemical equilibrium. We have already shown 
previously that the strangeness yield is following the pre- 
dicted QGP yield We also show that, when the en- 
ergy and baryon number deposited in the initial fireball 
drops below 20%, e.g., due to large impact parameter in- 
teractions, or/and small nuclei involved in the collision, 
the formation of the QGP phase becomes improbable. 

An important dynamical aspect of our experimental 
data analysis on which this work relies, is the sud- 
den disintegration (hadronization) of the fireball into the 
final state hadrons. This reaction mechanism is a pri- 
ori not very surprising, since a fireball formed in these 
collisions explodes, driven by internal compression pres- 
sure. However, we did not find in the particle production 
analysis as the particle source the expected chemically 
equilibrated, confined hadron phase. Furthermore, the 
particle production temperature (chemical freeze-out) we 
found corresponds to a deeply super-cooled state, which 
can be subject to mechanical instability |Q. We have 
checked elsewhere this assertion, applying the EoS de- 
veloped here in a detailed analysis of the fireball sudden 
break-up, and have found that it occurs at the condition 
of mechanical instability ||] . 

A consequence of the sudden fireball break up is that 
akin to the situation found in the study of properties 
of the early universe, particle abundances do not reach 
chemical equilibrium. Our analysis accounts in full for 
this important fact |||^]. For an appropriate criticism 
of the chemical equilibrium models, and a list of related 
work, we refer the reader to work of Biro || . 

In the next section, we address the thermal QCD inter- 
action coupling a s we will use. In section III, we define 
the equations of state and the parameters of our approach 
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and explore which is the best scheme for the extrapola- 
tion of the lattice data. In section IV , we present proper- 
ties of the QGP phase, relevant both to the study of the 
freeze-out conditions, and the study of the initial con- 
ditions reached in the collision. We present and discuss 
the comparison of the properties of the exploding fireball 
with those measured by means of hadron production in 
section [v| Our conclusions follow in section VI . 



II. QCD INTERACTIONS IN PLASMA 

The energy domain in which we explore diverse prop- 
erties of dense strongly interacting matter is barely above 
the scale 1 GeV and thus, in our consideration, an impor- 
tant input is the scale /x dependence of the QCD coupling 
constant, a s (/z), which we obtain solving 



da s 



-b Q al - ha 3 , + 



/3 2 port . 



(1) 



/3 2 ert is the beta-function of the renormalization group in 
two loop approximation, and 

11 - 2ra f /3 



b = 



2n 



h = 



51 - 19n f /3 
4^ 



/?2 GIt does not depend on the renormalization scheme, 
and solutions of Eq. (Q) differ from higher order renor- 
malization scheme dependent results by less than the er- 
ror introduced by the experimental uncertainty in the 
measured value of a s {u = M z ) = 0.118 + 0.001 - 0.0016. 
When solving Eq. (|]) with this initial condition, we cross 
several flavor mass thresholds and thus nf([i) is not a con- 
stant in the interval [i G (1,100) GeV. Any error made 
when not properly accounting for rif dependence on a 
accumulates in the solution of Eq. (|l|). In consequence, a 
popular approximate analytic solution shown as function 
of T, as dotted line in figure [0, 



dqL 



1 - 



2&i lni 



L ee ln( M 2 /A 2 ) , (2) 



with A = 0.15 GeV and rif = 3, is not precise enough 
compared to the exact 2-loop numerical solution. 

We show the numerical solution for a s (/i) in figure 
setting 



u = 2nT = kT/T c 



lGeV, 



where the solid line corresponds to a s (a = Mz) = 0.118, 
bounded by the experimental uncertainty Jl0| . We ob- 
serve that a s /n < 0.25 is relatively small. However, the 
expansion parameter of thermal QCD g = y/4:Tra s > 1 . 
Therefore even the small difference shown in figure [l] be- 
tween the approximation (dotted line) used in above ref- 
erences and the exact result is relevant. A large value of g 
has been a source of considerable concern about validity 
of perturbative expansion in thermal QCD O-Oj. 
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FIG. 1. a B (2ivT) for T c = 0.16 GeV. Dashed line: 
as(Mz) = 0.119; solid line = 0.118; dot-dashed line = 0.1156. 
Dotted line: approximate 2-loop solution given in Eq.(Q). 
Thin solid line: same as dotted, but extrapolating with rif = 5 
to a s (M z ) = 0.118. 

The analytical function, Eq. (||), is in fact a very good 
solution of Eq. (§). For A = 0.225 GeV and n f = 5, 
it agrees with numerical result for the entire range, 
M z = 92GeV> fi > M B ~ 4.5 GeV, in which n f = 
Const., here Mb is the bottom quark mass. However, 
the region of interest for the QGP equations of state we 
explore is 1 GeV< u < Mb, where nf varies and the ap- 
proximation, Eq. (H) is not adequate. The thin solid line, 
in figure ^, shows the behavior of the analytical solution 
with rif = 5 kept constant, and for A = 0.225 GeV which 
assures the boundary value a s (Mz) — 0.118. 



III. THE QUARK-GLUON LIQUID 

We now can define the quark-gluon liquid model which 
describes well the properties of QGP determined by 
lattice-QCD method. 

1. To relate the QCD scale to the temperature T = 1//3, 
we use for the scale the Matsubara frequency |15l : 



/i = 27r/r\/i 



~o In 2 



(3) 



This extension to finite chemical potential u qy or equiv- 
alently quark fugacity A q = exp// q /T, is motivated by 
the form of plasma frequency entering the computation 
of the vacuum polarization function Jl~6| ] - In principle, 
there should be in Eq. (||) also a contribution from the 
strange quark chemical fugacity A s expressed equiva- 
lently by the strange quark chemical potential u s . How- 
ever since strangeness conservation virtually assures that 
fi s ~ 0, or equivalently, A s ~ 1, we will not pursue this 
further. 

2. To reproduce the lattice results available at fi q = 
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pL we need to introduce, in the domain of freely mobile 
quarks and gluons, a finite vacuum energy density: 



B = 0.19 



GeV 



This also implies, by virtue of relativistic invariance, that 
there must be a (negative) associated pressure acting on 
the surface of this volume, aiming to reduce the size of the 
deconfined region. These two properties of the vacuum 
follow consistently from the vacuum partition function: 



lnZ v 



-BV0. 



(4) 



3. The partition function of the quark-gluon liquid com- 
prises interacting gluons, n q flavors of light quarks Jl7| , 
and the vacuum S-term. We incorporate further the 
strange quarks by assuming that their mass in effect re- 
duces their effective number n s < 1: 
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ci(^T) 4 
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where: 



ci 



15a s 



C-2 = l- 



47T 

50a s 
21tt 



(5) 



(6) 



c 3 = 1 

7T 



In figure the 'experimental' values are from numeri- 
cal lattice simulations of P/T 4 For practical reasons 
the lattice results for 'massless' 2 and 3-flavors were ob- 
tained with m/T = 0.4, which reduces the particle num- 
bers by 2%, and this effect is allowed for in the quark- 
gluon liquid lines (two flavors: dashed, three flavors: dot- 
ted) in figure ||. In the case 2+1 flavors, a renormalized 
strange quark mass m s /T = 1.7, the lattice input has 
been wP s /T = 1. This leads to a ~ 50% reduction in 
strange quark number, thus n s ~ 0.5 will be used in 
our study of QGP properties, assuming that strangeness 
is fully chemically equilibrated, unless otherwise noted. 
Thus, in general, we have rif = n q + n s ~ 2.5. 

The thin lines, in figure |[ correspond to the previously 
reported results for a quark-gluon gas fl4| , with first or- 
der QCD correction introduced using approximate value 
of a s (fi), Eq. (Q), without the vacuum pressure term. For 
T > 2T C , the disagreement is an artifact of the approxi- 
mation for a s . 

Lattice results were also informally reported Q for the 
energy density, with rif = 3, and are show in figure [jj, 
upper 'experimental' points. The energy density, 



CQGP 



dlnZ QGP (/3,A) 
Vd/3 



(7) 



is sensitive to the slope of the partition function. Wc 
show, in figure |, the energy density and pressure for 
two extreme theoretical approaches: the solid lines are 
for the model we described above (/i = 2irT, B = 
0.19 GeV/fm 3 ), the dotted lines are obtained in the ther- 
mal expansion, including all, up to fifth order, scale de- 
pendent terms obtained by Zhai and Kastening [ fl2[ , and 
choosing the scale fx = 2.6itT in free energy and in the 
coupling constant, so the result reproduces the pressure 
well (bottom dotted line in figure ||). 




Lattice-QCD result; 
compared with our quark-gluon liquid model (thick lines) and 
lowest order perturbative QCD using approximate a s , Eq. (^) 
(thin lines used in prior work dotted line 3 flavors, 

solid line 2+1 flavors, and dashed line 2 flavors. 

In the upper portion of figure |[ comparing these two 
theoretical approaches chosen to reproduce the pressure, 
we see a clear difference in the energy density. The 5th or- 
der energy density (dotted line) (g/^irf = (a s /n) 5 / 2 /32 
Jjjfl , disagrees with the lattice data also at high T, where 
these are most precise. Near to T ~ 1.5T C , the solid line is 
visibly better describing the pressure. In fact, systematic 
study of the behavior of the (g/4ir) n expansion reveals 
that higher order terms do not lead to a stable result for 
the range of temperatures of interest to us Jl^JlJ] . 

It is not uncommon to encounter in a perturbative ex- 
pansion a semi-convergent series. The issue then is how 
to establish a workable scheme. Our result implies that 
the choice of the Matsubara frequency 27rT as the scale 
(i of the running coupling constant has for yet unknown 
reasons, the effect to minimize the contribution of the 
sum of all higher order terms in the expansion Eq. (|TJ), 
even for moderate temperatures, once the vacuum pres- 
sure is introduced. The agreement between our model 
and the lattice calculations arises despite (or maybe be- 
cause) the fact that in our evaluation of the temperature 
dependence of the coupling constant a s (T) we ignored 
the dependence of (3^ ert on the ambient temperature scale 
jl8j . It seems that in fact the omission of this dependence 
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conspires with omission of the higher order thermal con- 
tributions in the partition function, yielding the remark- 
able agreement we see between the lattice results and our 
simple minded approach, which considers the first order 
corrections in a s only. 
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FIG. 3. Top: informal lattice-QCD results |J for e/T 4 at 
X q = 1, compared with our quark-gluon liquid model (solid 
line). Doted line is an alternative approach in which all terms 
in partition function are summed as given in and the 
scale is set at /i — 2.6irT. Bottom: published lattice pres- 
sure results j(| compared to the two approaches. Solid line: 
first order with the bag constant; dotted line: 5th order with 
/i = 2.6ttT. 

Given that our model adequately describes both pres- 
sure and energy density leads us to believe that we also 
have an appropriate extrapolation scheme for the QGP 
properties to finite baryon density - but we could not test 
if the coefficient C3 , Eq. (||) , describes well the behavior 
of the partition function for finite chemical potentials, as 
such lattice-QCD results are not available. 

In absolute terms, the model we now adopt for further 
study reproduces the lattice results well, at the level of a 
few percent. It can only be hoped that the lattice results 
have reached that level of precision. In our approach, the 
value of B we obtain and employ is entirely dependent 
on the quality of the lattice results. However, in another 
work la] , we have considered the dynamics of the expand- 
ing fireball and obtained an estimate for the value of B as 
defined here. We found that B > O.I7GeV/fm 3 , consis- 
tent with the value we employ here B — 0.19 GeV/fm . 



IV. PROPERTIES OF QGP-LIQUID 

We are now ready to explore the physical properties of 
the quark-gluon liquid. The energy density is obtained 
from Eq. (g) , recalling that the scale of the interaction is 
given by Eq. (§): 



e QGP = 4B + 3P QGP + A e + A q + A s 
A g =(b al + b ia l)^T\ 



A q = (boa^+hal) 
A s = (b a 2 s + b ia 3 s ) 



birn 



18 + 7T I ^ + 2^ 2 



7T V 2ir z 



18 



(8) 
(9) 

, (io) 

• (11) 



A convenient way to obtain entropy and baryon density 
uses the thermodynamic potential T: 

F{T,p,V) = _'L hlZ{ p ,X q ,V) QGP = -P QGP . (12) 
The entropy density is: 



S QGP 



dT 
~VdT ' 
(n q + n s )77r 2 



15 



c 2 T 3 + n q c 3 //qT + n s c 3 /i 2 T 



32tt' 



ciT 3 + {A g + A q + A s ) 



n 2 T 



tt 2 T 2 + /i 2 



(13) 



(14) 



Noting that baryon density is 1/3 of quark density, we 
have: 



Ph 



1 dT 



3 Vdfi q 



(15) 



+ -(A g + A q + A s ) 



Mq 



7T 2 T 2 + H 2 q 



(16) 



We show properties of the quark-gluon liquid in a wider- 
range of parameters in figure ||. We study the properties 
at fixed entropy per baryon S/b since an isolated ideal 
particle fireball would evolve at a fixed S/b. We consider 
the range S/b = 10 (at left for the top panel, baryo- 
chemical potential /Xb, and middle panel baryon density 
n/no, here no = 0.16/fm 3 , and bottom left for the energy 
per baryon E/b) in step of 5 units, up to maximum of 
S/b = 60. The highlighted curve, in figure |], is for the 
value S/b = 42.5, which value follows from earlier study 
of hadronic particle spectra |||. The dotted line, at the 
minimum of E/b\s/bi is where the vacuum and quark- 
gluon gas pressure balance. This is the equilibrium point 
and indeed the energy per baryon does have a relative 
minimum there. 

Unlike the case for an ideal quark-gluon gas, the lines 
of fixed S/b, seen in the top panel of figure [| are not corre- 
sponding to /J,b/T = Const., though for large T and small 
Hb they do show this asymptotic behavior. Since little 
entropy is produced during the evolution of the QGP fire- 
ball, the thick line in the lower panel of figure ^ describes 
the approximate trajectory in time of the fireball made in 
Pb-Pb interactions at the projectile energy 158A GeV. 
We stress that there has been no adjustment made in any 
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parameter to bring the earlier determined 'experimental' 
points shown in figure |] into the remarkable agreement 
with the properties of equations of state of the quark- 
gluon plasma, also seen in figure pL 
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FIG. 4. From top to bottom: fib, n/no and E/b; lines 
shown correspond to fixed entropy per baryon S/b — 10 to 
60 by step of 5 (left to right). Thick solid lines: result for 
S/b = 42.5. Limits: energy density £ q , g = 0.5GeV/fm 3 
and baryo-chemical potential fit = 1 GeV. The experimental 
points denote chemical freeze-out analysis result || , discussed 
in section ^. 

The trajectories at fixed energy per baryon E/b are 
shown in the T~X q plane in figure [^, for the values 
(beginning at right) E/b = 2.5 to 9.5 GeV by step 
of 1. The highlighted curve corresponds to the value 
E/b = 7.8 GeV which is the local intrinsic energy con- 
tent of the hadronizing QGP fireball formed at SPS Pb- 
Pb interactions at the projectile energy 158A GeV [||. 
Dotted line, in figure |[ corresponds to P = 0, the solid 
line that follows it is the phase transition line where the 
QGP pressure is balanced by pressure of the hadron gas. 



This is the condition at which the equilibrium transition 
would occur in a slowly evolving system, such as would 
be the early universe. The properties of hadronic gas 
are obtained resuming numerically the contribution of 
all known hadronic particles including resonances, which 
effectively accounts for the presence of interactions in the 
confined hadron gas phase p9[. 
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FIG. 5. Contours of constant energy per baryon in QGP 
in the T-\ plane: From right to left E/b = 2.5 to 9.5 GeV 
by step of 1, thick solid line is for E/b = 7.8 GeV. Dotted line 
corresponds to P — 0, above this, the solid line is the phase 
transition where the QGP pressure is balanced by pressure 
of the hadron gas. The experimental point denotes chemical 
freeze-out analysis result |3|, see section |v| Bottom dashed 
line boundary is at energy density 0.5GeV/fm 3 . 

In relativistic nuclear collisions, the formation of equi- 
librium state competes with the evolution of the fireball. 
The slowest of all the equilibria is certainly the chem- 
ical equilibration of strange quarks ^0|. Next slowest 
is the equilibration of light quarks. The chemical relax- 
ation time constant for the production of light quarks 
has been obtained in the first consideration of chemi- 
cal QGP equilibration, see figure 2 in pfj| : Tca^qq{T — 
250 MeV) = 0.3 fm. The chemical equilibration of glu- 
ons has also been questioned, and is found to be slow, 
when only 2G — > 3G processes are allowed |2l| . Since 
gluon fusion processes are proportional to the square of 
gluon abundance, Tca^qq{T — 250 MeV) — > 1.2 fm for 
processes occurring when gluons are at 50% of equilib- 
rium abundance. This explains why quarks trail gluons 
in the approach to equilibrium, which are struggling to 
equilibrate by multi-gluon production processes p^| . It 
is generally assumed that the approximate thermal (ki- 
netic) equilibrium is established much faster. The mech- 
anisms of this process remain under investigation p3| ] . 

In figure ^, a case study how the chemical equilibration 
of quarks cools the gluon-chemically equilibrated fireball 
is presented. The change of both the initial temperature 
Ti (upper panel) and X q (lower panel) as function of n/, 
for E/b = 9.3 GeV, the total final energy content of Pb - 
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Pb collisions at CERN, is presented. We see that in the 
initial state an equilibrated glue phase at T ~ 270 MeV 
may have been reached. Full chemical equilibrium would 
correspond to T > 230 MeV, but more likely this value is 
somewhat reduced due to flow dilution that has occurred 
in the process of chemical equilibration. 

300 I i i i i | i i i i | i i i i | i i i i | i i i i _ 




FIG. 6. The initial conditions as function of the number 
of Fermi degrees of freedom nf. Upper panel for temperature 
and lower for quark fugacity, at E/b — 9.3 GeV. 

One of the interesting quantities is the QGP energy 
density which we show in figure @, for both fixed E/b 
(top) and fixed S/b bottom, for fully equilibrated con- 
dition n q = 2,n s = 0.5. We see that for E/b > 6 GeV 
and S/b > 25 the influence of finite baryo-chemical po- 
tential is minimal and the lines coalesce. In other words, 
at conditions we encounter at SPS, we can correlate the 
energy density with temperature alone e ~ aT 4 , as seen 
for nf = 3 > 2.5 in figure ||. In figure 0, we see that 
the energy density 3GeV/fm 3 is established when the 
temperature in the equilibrated fireball equals 212 MeV. 
Considering the observed high inverse slopes of strange 
particles, one can assume that the the plasma phase, 
before it reached full chemical equilibrium, has been at 
about T ~ 250 MeV and the energy density in this state 
has most likely been still higher, above 4GeV/fm 3 . We 
believe that our evaluation of the properties of the QGP 
liquid is, at this range of temperature T ~ 1.5T C , reliable. 

As a final step in the study of the properties of the 
QGP liquid, we consider the conditions relevant for the 
formation of the QGP, and consider the behavior for 
nf = 1. We show, in figure ||, lines of fixed energy per 
baryon E/b =3, 4, 5, 6, 8, 10, 20, 50 and 100 GeV, akin to 
results we have shown for n/ = 2 + 0.5, in figure |[ The 
horizontal solid line is where the equilibrated hadronic 
gas phase has the same pressure as QGP-liquid with semi- 
equilibrated quark abundance. The free energy of the 
QGP liquid must be lower (pressure higher) in order for 



hadrons to dissolve into the plasma phase. The dotted 
lines in figure ||, from bottom to top, show where the 
pressure of the semi-equilibrated QGP phase is equal to 
T) = 20%, 40%, 60%, 80% and 100%, t] being the 'stop- 
ping' fraction of the dynamical collisional pressure p4| : 




I i i i i i i i i i I i i i i i i i i i I 

0.1 0.2 0.3 

T [GeV] 

FIG. 7. Energy density in QGP as function of tempera- 
ture. Top panel: for a fixed Energy per baryon E/b = 2 to 
15 GeV by step of 1; bottom panel for fixed entropy per baryon 
S/b = 10 to 60 by step of 5. Boundaries are e = 0.5 GeV/fm 3 
at the bottom and yn, = 1 GeV at the top. 

The rationale to study, in figure pi lines at fixed E/b 
is that, during the nuclear collision which lasts about 
2i?7v/7i2c ~ 13/18 fm/c, where jl is the Lorentz factor- 
between the lab and CM frame and Rn is the nuclear 
radius, parton collisions lead to a partial (assumed here 
to be 1/2) chemical equilibration of the hadron matter. 
At that time, the pressure exercised corresponds to col- 
lisional pressure P co i . This stopping fraction, seen in 
the transverse energy produced, is about 40% for S-S 
collisions at 200,4 GeV and 60% for Pb-Pb collisions at 
158 A GeV. If the momentum-energy and baryon number 
stopping are similar, as we see in the experimental data, 
then the SPS collisions at 160-200A GeV are found in the 
highlighted area left of center of the figure. In the middle 
of upper boundary of this area, we would expect the be- 
ginning evolution of the thermal but not yet chemically 
equilibrated Pb-Pb fireball, and in the lower left corner of 
the S-S fireball. We note that the temperature reached 
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in S-S case is seen to be about 25 MeV lower than in 
the Pb-Pb case. The lowest dotted line (20% stopping) 
nearly coincides with the non-equilibrium phase bound- 
ary (solid horizontal line, in figure |^) and thus we con- 
clude that this is, for the condition rif = 1, the lowest 
stopping that can lead to formation of a deconfined QGP 
phase. Such a low stopping would be encountered possi- 
bly in lighter than S-S collision systems or/and at large 
impact parameter interactions of larger nuclei. 

The highlighted area, right of center of the figure [|, 
corresponds to the expected conditions in Pb-Pb colli- 
sions at 40A GeV. If we assume that the stopping here 
is near 80%, then the initial conditions for fireball evo- 
lution would be found towards the upper right corner of 
this highlighted area. We recognize that the higher stop- 
ping nearly completely compensates the effect of reduced 
available energy in the collision and indeed, we expect 
that we form QGP also at these collision energies. It is 
important to realize that we are entering a domain of 
parameters, in particular X q , for which the extrapolation 
of the lattice results is not necessarily reliable, and thus 
our equations of state have some systematic uncertainty. 
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FIG. 8. Contours of energy per baryon in QGP in the 
T-\ q plane for n/ = 1: From right to left E/b=3, 4, 5, 
6, 8, 10, 20, 50 and 100 GeV . Thin, nearly horizontal line: 
hadronic gas phase has the same pressure as the QGP-liquid 
with semi-equilibrated quark flavor. Dotted lines from bottom 
to top: pressure in QGP liquid equals 20%, 40%, 60%, 80%, 
and 100%, of the dynamical collisional pressure. 



V. EXPERIMENTAL DATA ANALYSIS AND 
QGP-EOS 

A full account of our prior analysis of the 158 A GeV 
Pb-Pb collision system has appeared HQ. We briefly 
summarize the results that we require for the study of 
QGP properties at fireball breakup. In table |[ in up- 
per section, we present the parameters Tf (the chemical 
freeze-out temperature), v c (the collective flow velocity 



at sudden breakup), X q (the light quark fugacity), A s 
(the strange quark fugacity), j q (the light quark phase 
space occupancy) , 7 S (the strange quark phase space oc- 
cupancy) . These are derived from analysis of all hadrons 
excluding f2 and Q, which data points are not following 
the same systematic production pattern. These param- 
eters characterize completely the physical properties of 
the produced hadrons, and these properties are shown in 
the bottom section of table |. 

In the heading of the table, the total error, x 2 is shown, 
along with the number of data points N, parameters p 
and data point constraints r. The confidence level that is 
reached in our description is near or above 90%, depend- 
ing on scenario considered. The scenarios we consider are 
seen in the columns of table [|: an unconstrained descrip- 
tion of all data in the first column, constraint to exact 
strangeness conservation in the observed hadrons, second 
column. Since in both cases the parameter j q assumes 
value that maximizes the entropy and energy content in 
the pion gas, we assume this value in the so constrained 
third column. 

We can now check the consistency between the statis- 
tical parameters (top panel of table |) and the physical 
properties of the fireball (bottom panel of table |) which 
are maintained in the process of hadronization. We note 
that the energy shown in this table, is the intrinsic energy 
in the flowing frame. The CM-laboratory energy includes 
the kinetic energy of the flow and thus is greater, to be 
obtained multiplying the result shown in table || by the 
Lorentz factor 7 = — v% = 1.19. Thus the initial 

value of the energy per baryon that the system has had 
before expansion started has been E°/b ~ 9.3 GeV, as 
used in figure ^| and in the estimate presented in figure ||. 

TABLE I. Results of study of Pb-Pb hadron production 
||: in the heading: the total quadratic relative error Xti 
number of data points N, parameters p and redundancies 
r; in the upper section: statistical model parameters which 
best describe the experimental results for Pb-Pb data. Bot- 
tom section: specific energy, entropy, anti-strangeness, net 
strangeness of the full hadron phase space characterized by 
these statistical parameters. In column one, all statistical 
parameters and the flow vary. In column two, we fix A s by 
requirement of strangeness conservation, and in column three, 
we fix 7 9 at the pion condensation point "f q — 





Pb|„ 


Pb|; b 


Pb|J c 


Xt! N;p;r 


2.5; 12; 6; 2 


3.2; 12; 5; 2 


2.6; 12; 5; 2 


T t [MeV] 


142 ± 3 


144 ± 2 


142 ± 2 


Vc 


0.54 ± 0.04 


0.54 ± 0.025 


0.54 ± 0.025 




1.61 ± 0.02 


1.605 ± 0.025 


1.615 ± 0.025 


A, 


1.09 ± 0.02 


1.10* 


1.09 ± 0.02 




1.7 ± 0.5 


1.8± 0.2 


7 c* = e m*/w> 


7s/7 9 


0.79 ± 0.05 


0.80 ± 0.05 


0.79 ± 0.05 


Ef/B 


7.8 ± 0.5 


7.7 ± 0.5 


7.8 ± 0.5 


Sf/B 


42 ± 3 


41 ± 3 


43 ± 3 


s f /B 


0.69 ± 0.04 


0.67 ± 0.05 


0.70 ± 0.05 


(s f -s f )/B 


0.03 ± 0.04 


0* 


0.04 ± 0.05 
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In the bottom panel in figure kl we saw that the 
Temperature Tt = 143 ± 3 MeV and energy per baryon 
E/b = 7.8 GeV where just at S/b = 42.5 seen table @. 
Similarly, in the top panel, the baryo-chemical potential 
fi b = 3Tf In X q = 204 ± 10 MeV is as required for the 
consistency of QGP properties. A similarly embarrass- 
ing agreement of the hadron yield analysis results with 
properties of the QGP fireball is seen in figure |, but 
that is just a different representation (at fixed E/b) of 
the result we saw at fixed S/b, in figure ||. However, 
the importance of this result is that the plasma breakup 
point appears well below the phase transition tempera- 
ture line (thin horizontal solid line). As this discussion 
shows, the properties of the QGP liquid at the point of 
hadronization are the same as found studying the prop- 
erties of hadrons emerging from the exploding fireball. 
Both the specific energy and entropy content of the fire- 
ball are consistent with the statistical parameters Tf and 
Hb according to our equations of state of the quark-gluon 
liquid. The freeze-out point at T ~ 143 MeV, seen in bot- 
tom panel of figure ||, corresponds to an energy density 
£f — 0.6 GeV/fm 3 . This is the value for the super-cooled 
plasma, the equilibrium phase transition occurs at twice 
this value, e p ~ 1.3GeV/fm 3 . 

We can also evaluate the hadronic phase space en- 
ergy density. We need to introduce the excluded vol- 
ume correction |2j|. Considering that the point parti- 
cle phase space energy density e pt — l.lGeV/fm 3 , we 
obtain ehg — 0.4GeV/fm 3 , using the value of B = 
0.19GeV/fm 3 . Taking into account the numerous un- 
certainties in the understanding of the excluded volume 
effect, we conclude that the fireball energy density is com- 
parable in magnitude to the energy density present in 
hadronic phase space. We also note that the equilibrium 
phase transition curve we had presented in figure ^| was 
computed without the excluded volume effect. When al- 
lowing for this, the solid line (critical curve) moves about 
half way towards the dotted line where -Pqgp = 0. 

For a vanishing baryo-chemical potential, X q = 1, we 
determine in figure ^, that the phase transition temper- 
ature for point hadrons is T p ~ 172 MeV. The super- 
cooled P = temperature is at T c — 157.5 MeV (dotted 
line at X q = 1, in figure |^), and an expanding fireball 
can super-cool to as low as T ~ 145 MeV, where the me- 
chanical instability occurs We also have seen this 
result in the bottom panel of figure [| where the doted 
line, corresponding to pressure zero, is to the right of the 
'experimental' point. 

The collective velocity v c of the exploding matter 
remains large in deeply super-cooled conditions even 
though the expansion slows down and kinetic energy is 
transfered back from flow to thermal component, once 
the dotted P = line is crossed, as can be seen in fig- 
ure ^. Interestingly, our central value v c = 0.54 is the 
velocity of sound of the exploding fireball: 

v ° = £ • (i7) 

° £ S/b 



The theoretical line shown in figure |9j is computed for 
S/b = 42. 




0.5 

100 200 300 400 500 

T [MeV] 

FIG. 9. The velocity of sound of quark-gluon liquid at 
S/b — 42 . The dotted line corresponds to the value of the 
sound velocity of an ideal relativistic gas, v s = l/Vo. 

There is yet more evidence for the QGP nature of the 
fireball. It has been argued before that the values of the 
three other chemical parameters A s , 7 9 and 7 S suggest as 
source a deconfined QGP fireball source ||: 

a) The value of strange quark fugacity A s can be ob- 
tained from the requirement that strangeness balances, 
(N s — N s ) — , which for a source in which all s an 
s quarks are unbound and thus have symmetric phase 
space, implies A s = 1 . However, the Coulomb distortion 
of the strange quark phase space plays an important role 
in the understanding of this constraint for Pb-Pb colli- 
sions, leading to the Coulomb-deformed value A s ~ 1.1, 
which is identical to the value obtained from experimen- 
tal data analysis A" = 1.09 ± 0.02 . 

b) The phase space occupancy of light quarks j q is, 
before gluon fragmentation, near or at the equilibrium 
value 7 g = 1. However, as measured by hadron abun- 
dances it is expected to significantly exceed unity to 
accommodate the contribution from gluon fragmenta- 
tion into light quark pairs. There is an upper limit: 
7<j < lq = e m,r / 2T/ ~ 1.67, which arises to maximize 
the entropy density in the confined hadron phase. 

c) The strange quark phase space occupancy 7 S can be 
computed within the framework of kinetic theory and 
is mainly influenced by strangeness pair production in 
gluon fusion [pb[ , in early stages of the collision at high 
temperature, and by dilution effect in which the already 
produced strangeness over saturates the 'thiner' low tem- 
perature phase space. Moreover, some gluon fragmenta- 
tion also enriches 7 S as measured by hadron abundance. 
We note that some earlier studies implicitly address the 
parameter j s /jq which therefore is stated in table Q. 
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It is worthwhile to recall here that the strangeness 
yield, s/b — s/b ~ 0.7, predicted early on as the re- 
sult of QGP formation [^0| is also one of the results of 
data analysis seen in table |. This result is found with 
modern kinetic studies of strangeness production in QGP 

IS- 

VI. DISCUSSION AND CONCLUSIONS 

We have presented properties of the quark-gluon 
plasma equations of state, and have made several appli- 
cations pertinent to the physics of relativistic heavy ion 
collisions. These included a study of the initial attainable 
conditions reached in Pb-Pb collisions at 158AGeV pro- 
jectile energy, as well as an evaluation of the hypothesis 
that direct hadron emission occurs from a disintegrating 
QGP fireball. 

We have based our description of the QCD matter on 
a form obtained in perturbative expansion of thermal 
QCD, wherein we have introduced a nonperturbative re- 
summation of the thermal interaction strength. We have 
also introduced the non-perturbative vacuum properties. 
The two parameters of our model are the vacuum pres- 
sure B and the scale //(T) of the QCD coupling strength. 
These were chosen to reproduce the latest lattice results 
for pressure and energy density obtained at zero baryon 
density. For study of finite baryon density we rely on 
the perturbative behavior of the interacting quark-gluon 
gas. This is developed within the same computational 
approach in which the pressure and energy density are 
obtained correctly. Fortuitously, it turns out that the 
QCD interaction effects are smallest for the baryon den- 
sity, and thus, a priori, are most reliably described. Even 
so, the magnitude of these effects increases to the level 
of 45% at the phase boundary. 

A serious limitation of our approach is that we rely 
completely on the current lattice results. There are open 
questions about the precision of lattice results we model. 
The requirement for continuum extrapolation Q , and the 
use of the relatively large m/T — 0.4 make our results 
uncertain at the level of 10%. Moreover, in the deeply 
supercooled region we do not have lattice results available 
and the behavior of QGP we consider arises solely from 
the analytical properties of our QGP model. 

In the region of statistical parameters of interest in 
the discussion of the SPS experimental data, the detail 
how the baryo-chemical potential enters the interaction 
scale /i of the strong interaction does not matter since 
for {nT c ) 2 ~ 50fi q the ^-contribution to the interac- 
tion scale is negligible compared to the T-contribution. 
This, however, means that our successful description of 
the SPS-experimental results does not imply that our 
model of QCD matter can safely be used in the study 
of properties of quark star matter, or the fireball matter 
made at AGS energies, where T c ~ p, q [p6[ . 

The validity of our equation of state model is better 



established in the temperature range pertinent to study 
of the initial state of the QGP fireball formed at SPS. 
We have found that during the pre-chemical equilibrium 
stage of light quarks, the so called hot glue scenario 
p2|, the QGP plasma at SPS has been formed at about 
T ~ 250-270 MeV, see figures | and § By the time 
light quarks have also chemically equilibrated, the tem- 
perature decreases to just above T ~ 210 MeV, and the 
energy density in such a QGP fireball is at 3GeV/fm 3 , 
as can be seen in the lower panel of figure 0. 

We have studied the requirement for the minimal col- 
lision energy, which permits the formation of deconfincd 
fireball as function of 'stopping' in the collision. We have 
shown that a significant change in the reaction mecha- 
nism of colliding nuclei can be expected when the opac- 
ity is at or below 20%, see figure ||. At that point, for- 
mation of a deconfined phase fireball is not assured as 
is seen comparing in the T-\ q plane the boundary be- 
tween hadron phase (solid horizontal line in figure |SJ) and 
chemically non-equilibrated QGP initial conditions (low- 
est dashed horizontal line in figure ||) . Such a low stop- 
ping is expected for sufficiently small collision systems, 
with participant mass below that seen in central S-S in- 
teractions. The results shown in figure || imply that the 
QGP-liquid phase should be formed for Pb-Pb collisions 
at all collision energies accessible to SPS. In fact, were 
it not for the uncertainties inherent in the extrapolation 
of lattice results to high baryon density reached at AGS, 
we could argue that deconfinement also occurs for Au-Au 
collisions in the high AGS energy range |^7| . 

We have also compared the properties of the QGP 
phase to conditions present at time of hadronization, 
as obtained within an analysis of hadron production |jj . 
This analysis determines in a first step a set of statisti- 
cal parameters which describe well experimental hadron 
multiplicity results. Since this set of parameters also 
characterizes the phase space of all hadronic particles, 
in a second step one can estimate the energy, baryon 
number and entropy content contained in all hadrons 
produced. A comparison to the QGP equations of state 
can be made, assuming that in the hadronization of the 
deconfined matter there has been no reheating, and no 
shift in chemical equilibrium properties. This is the sud- 
den hadronization. In this case we can compare final 
state fireball properties with the behavior of QGP mat- 
ter evaluated at same values of T, X q , X s . 

However, the momentum distribution of final state 
hadrons could be subject to modification after QGP fire- 
ball breakup, which is usually expressed by introducing 
two sets of statistical parameters, the chemical freeze-out 
(particle abundance freeze-out) and thermal freeze-out 
(spectral shape) temperatures and chemical potentials. 
In our comparison of QGP properties with fireball prop- 
erties we rely on results obtained from an analysis of 
particle multiplicity ratios, which reduce the impact of 
post-hadronization system evolution. 

At the time of hadronization of a QGP fireball gluons 
contribute to form a chemical nonequilibrium excess of 
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hadrons. Thus even if the hadronizing QGP is (near) 
chemical equilibrium, the confined phase in general will 
require allowance for chemical nonequilibrium. A non- 
equilibrium analysis is more general than the equilibrium 
models |l, and always describes the hadron production 
data better ||. This result in itself constitutes evidence 
for the presence of primordial deconfmcd phase. Here, 
we have for the first time demonstrated that there is a 
good agreement between the energy, entropy and baryon 
number content of the collision fireball and the properties 
expected from the study of supercooled QGP phase. 

Although this comparison has produced a very good 
agreement between the fireball properties and the QGP 
properties, this could be just a coincidence. Firstly, we 
do not fully understand the limitations in the description 
of the fireball matter. The hadron phase space is used 
with statistical populations allowing for presence of many 
very heavy hadron resonances. This hypothesis has not 
been experimentally verified at this time. In fact it is to 
be expected that a more complex pattern of formation 
of hadronic particles arises, and we estimate the system- 
atic error in the computed physical properties of hadron 
phase space to be at least 15%, arising from the assumed 
populations of unobserved hadronic states. Because of 
the nature of the QGP equations of state the 'measure- 
ment' point in figure ^ presented at the values of fib and 
T found from hadron multiplicity analysis would be com- 
patible with a range 35 < S/b < 50. While it is nice to 
see the evaluated hadron abundance property to be just 
it in the middle of this QGP range, we must keep in mind 
that a wide range of entropy values is permissible. 

It must be clearly said that the fireball properties we 
study are incompatible with properties of chemically and 
thermally equilibrated confined hadron gas. Perhaps the 
simplest way to see this is to realize that applying such 
a model to describe strange hadrons gives a set of pa- 
rameters which fail to describe the total hadron yield by 
many (10-15) standard deviations. The key point of our 
study is that we have shown that a natural way to ex- 
plain consistently all hadronic production data obtained 
in 158^4 GeV Pb-Pb collisions is to invoke as reaction 
picture the formation, and sequel sudden hadronization, 
of a QGP fireball. 
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